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Abstract 

In this paper wc first obtain local contraction results in the H m x 
ff 2m -norm with m > 1 with respect to time and space for a local 
scheme. Especially, we show that the local scheme preserves C m (~1 H m 
regularity together with polynomial decay of order m of the data at 
each time step for each m > 2. For an extended controlled scheme 
with local limit functions / — > v[' p ' ( (Z,.) = v?' (I, .) + r\ we observe 
linear growth with respect to the time step number I and with respect 
to the H m x H 2m -norm. We simplify the controlled scheme considered 
in pQ and [2J. For the simplified controlled scheme we observe linear 

growth of the squared local limit functions I — > (v?' l (l, .)j measured 

in H m x H 2m norm . This leads to a global linear bound of the Leray 
projection term. Furthermore, the schemes discussed in [1] and [2] 
are simplified in the sense that the estimates are achieved without the 
use of some properties concerning the adjoint of a local fundamental 
solutions with variable drift terms. We note that the pointwise and 
absolute convergence of the local functional series and their first order 
time derivatives and their spatial derivatives leads to a constructive 
approach of local classical solutions. 

2000 Mathematics Subject Classification. 35K40, 35Q30. 

1 Introduction 

First, we recall a naive simple scheme in brief which leads to loca solutions. 
For time step I > 1 and small < pi ~ -j- we consider a time-local functional 
scheme for functions v?' l,k ' , 1 < i < n.k > with limit 



V; 



?' l :=ve' l -\i-i,.) + YM> l > k , (1) 



k=l 



where for 1 < i < n we have v?' 1 1 (l — 1, .) € H 2 n C 2 , and Sv p,l ' k = 
v P' l ' k — y P' l ' k ~ l are functions denned on the domain [I — 1, 1] x M. n along with 
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v P,l, i .— y P' 1 1 ^ ana > w ]2 ere the functions v P,l ' k satisfy the equations 

( dvf' k d 2 vf' k , on p,*,fc-l ^f' fc 

-fc^-pf / Ei=i-^3-+PiEi= 1 «i" -*r = 

« Elm=i k« {£rK n (x - yj) E], m=1 ( § ^ 9 -i^ 1 ) (r,y)dy, 
VP' l ' k (l-l,.) =\P' l - 1 (l-l,.) 

(2) 

on the domain [I — 1, 1] x M n for I > 1. Here, iT n is the fundamental solution 
of the Laplacian of dimension n. We are interested in the case n > 3. For 
/ = 1 we define v p '' _1 (Z — 1, .) = h(.). In this part III of articles on the mul- 
tivariate Burgers equation and the incompressible Navier Stokes equation 
we make four contributions. First we define a simplified scheme which will 
allow us to simplify the proof of local contraction in .ff 2 -based spaces in the 
sense that the argument based on the adjoint of the fundamental solution is 
avoided. Second, we improve the result concerning the local contraction by 
showing that a certain kind of polynomial decay is preserved together with 
higher order regularity. Third, we explicitly show why the scheme proposed 
is an approach of constructing classical solution rather than a mere numer- 
ical scheme. And fourth, we show that an extension of the scheme with a 
simplified control function leads H m x H 2m -base& linear bound of the growth 
of the solution with respect to the time step number I > 1 holds also for the 
squared norm of the controlled solution function, such that the controlled 
scheme is really global. The price to pay for a simplified control function is 
that we have only a linear bound and a linearly decreasing time-step size. 
We mentioned earlier in pQ that a controlled scheme of the incompressible 
Navier Stokes equation allows for a constant step size and for a numerical 
stabilization of computations. Here we propose an alternative proof method 
for the sake of global existence and regularity. However, we shall see that 
the essential step is the local contraction estimate. 



2 A simplified scheme 

For m > 2 and time step I > 1 and small < p\ ~ \ we consider a time-local 
functional scheme 

oo 

v *, P ,i :=v *, P ,i-i + j2sv*' p ' l > k , (3) 

k=l 

where for 1 < i < n v*'?' 1 ' 1 ^ — 1, .) G H m n C m , and 5v*<P< l < k = vf' k - 
v *,p,l, -l a iong with v* ,p,l ~ l := v*' p ' l ~ l , and where the functions v*' p ' l ' k sat- 
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isfy the equations 



f d 2 v*' p > l ' k _ s-^n *,p,l,k-ldv*> p ^ k - 1 
— Pl U l^j=l dxj ~ ~Pl A?'=l V j Wj 

( +pi El m =i / R - (irM* ~ v)) El m =i ("'v """C ) (r,y)dy, 

v *W,fe(/ - 1,.) = v*'^- 1 ^ - 1,.). 

(4) 

We are interested in the case n > 3, and for I = 1 we define v*^' _1 (Z-l, .) = 
h(.). Especially, we are interested in the case n = 3 of course. If S demotes 
the right side source term of the first equation of (HJ), then we see that we 
have the representation 

v*' p ' 1 '\t,x) = ; R8 v*'P' l -\l - l^G^x- y)dy 

+ Si-i Imp s '( s >y) G! K* - s,x - y)dyds 

where Gi is the fundamental solution of the heat equation ^ — pii/Au = 
on [/ — 1,/] x IR n . Compared to the previous scheme we observe that the 
latter expression is a convolution. A priori estimates are easier at hand, 
and we do not need the adjoint of the fundamental solution in order to shift 
derivatives. Nevertheless, keeping book of the additional source term we 
may adopt arguments of part I and part II and arrive at the same result by 
simpler considerations. 

3 Statement of local contraction result (strong form) 

The functions v? ,k ' 1 are define locally on [I— 1, 1] x M n . For each I > 1 we may 
consider the trivial extension (denoted by the same symbol 

V P, ,1 Q £ no ^ a tion for the sake of simplicity) which is defined to be zero on 
the complementary domain (R \ [I — 1, 1]) x M n x R. We may measure these 
extensions in I? (M n+1 ) = L 2 (R) x I? (W n ) obviously, where the latter set 
product notation indicates that time and space may be measured naturally 
up to different orders of regularity. We have the following local contraction 
result. 

Theorem 3.1. If v?'~ l (l — 1, .) G H 2 for 1 < i < n, then for small pi > 
with pi ~ j in the scheme above we have 

max ie{li ... in} \5vf' k \ H i xH 2 

(6) 

< imax i6{1] ... jn} \5vP' l ' k 1 \ H i xH 2. 
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Furthermore we may choose pi above such that in addition we have 
maxjgji,.. ,„} \5v?' 1,1 (t, .)\mxW 

(7) 

= max ie{li ...,„ } \vP' l '° - vP' 1 1 \ H i xH m < 

We shall see that this is indeed essential together with linear growth 
in order to get a global scheme, since for the coefficient function evaluated 
at t > the embedding C a C H 2 for dimension n = 3 (uniformly with 
respect to time t > 0) ensures that given the solution this same solution can 
be represented in terms of fundamental solution of certain scalar parabolic 
equations which involve the solution in the first order terms. However, in 
this paper we shall see that we can even improve this as follows. We say 
that a function / : R n — > R with / € C m (i.e., / has continuous partial 
derivatives up to order m) has polynomial decay of order m > 1 up to the 
derivatives of order m > 1 if for all multivariate partial derivatives D%f with 
\a\ = XT=i \ a i\ an d \ a \ < m we have for all x S R n 

\D2f(x)\ < - Cq , (8) 
1 x y ' l ~ 1 + \x\ m y 1 

for some C a < oo. Now we can formulate our improvement of theorem 13.11 

Theorem 3.2. If for 1 < i < n the function v* ,p,l ~ l (l — 1,.) G C m has 
polynomial decay of order m > 1 up to the derivatives of order m > 1, 
then for all k > we have that v*' p,k,l (T, .) for r G [I — 1,1] is in C m and 
has polynomial decay of order m > 1 up to the derivatives of order m > 1. 
Furthermore for small j < pi > in the scheme we have 

r *,p,i,ki 

max ie{l,-,n} \° v i \H™xH 2 ™ 

(9) 

,1 I j *,p,l,k— li 

and «;e may choose pi above such that in addition we have 
max ie{i,- ,n} l^ l *' P " 1 |j3" m xi? 2m 

I *,p,i.O/ \ *,P,l— 1/7 1 M ^ \ 

= sup T e[l-i,q niax ie{i,-,n} l«i ( r > •) ~ V (* ~ !• Oli^xfP™ < jyf- 

(10) 

4 Proof of theorem 13.21 

In order to prove the contraction property ([6]) we first consider representa- 
tions of 5v?' ' = vf ' ' — v?' ' 1 in terms of the fundamental solution Gi 
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defined above. First observe that the equation for 5vf ,l ' k is 

' d 5v*' p ' l ' h ™ d 2 5v*' p ' l ' h 

En *,p,l,k-ld6v*' p ' l ' k ~ 1 ^ n <. *,p,l,k—l dv*' p ' l ' k ~ 1 

j= i Vj — ^ Pi E j= i Svj — 

-A^it (" r: C """C ' ) (T.y)^, 

Sv*>P> l > k (l -!,.) = 0. 

(11) 

Again this equation is considered on the domain [I — 1,1] x R n . Next in 
terms of the fundamental solution Gi, and for all 1 < i < n we have the 
representation 

8v* ,p,l ' k (T,x) = 

+Plll-ilnn (-E"=i^*' Mfc_1 ^ i i- — ) (s,y)Gi{T-s;x-y)dyds 

+Pj/i-i/k» (" Ei=i«i' P ' , ' fc ~ 1 ^K; ) (s,y)Gi(T-s;x-y)dyds 

+PI Ii-i I«n J R n (£-K n (z ~ y)) x 

/ 8t)v*' p ' l ' k ~ 1 /a *,p,!,fe-l a *,p,i,fe-2\\ 

x (E^=i (^fe— + Hi-)) (*, y)G,(r - s - *)4ycfech, 

(12) 
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where for convenience we rewrite the difference of the Leray projection terms 
observing that 



, a *, P ,l,k-l Qy*^' 1 ^- 1 



2—/j,m=l 



dxj 



dx r , 



En 
j,m=l 



dv 



dxi 



dx n 



dx , 



8x r , 



En 
j,m=l 



dvt 



p,l,k-2 Q v *' p ' l ' k ~ 2 



dXn 



dx„ 



En 
7,m=l 



dSv 



*, P ,l,k-i dv*' p ' l ' k - 1 



dxj 



8x r , 



+ V" 
1 <j,m=l 



dvt 



P ,l,k-2 dSv*'"' 1 ^' 1 ' 



dxj 



dx r , 



(13) 



En 
rr, 



m,j=l 



"" v i dv 

3 uu rr, 



p,l,k-l 



dx,. 



dxj 



En 
3, 



771=1 



B *,p,l,k-2 d5v. 



,p,l,k-l ' 



dxj 



dx r , 



sr^n 

2-^777,7 = 1 



d8v 



.,p,!,fc-l 



dx„ 



dv 



:, P ,l,k-l 



dxj 



,p,l,k—2 



4.1 L 2 x L 2 -estimates 

We shall use Young's inequality. Let g G L 1 ([a, a + 1] xR n ) and / G 
L 2 ([a,t] x R n ) for some + + l>t>a>0. Basically in our the Young 
inequality is the observation that 



ll /m« /( s > y)9(t ~s,x- y)dydi 



\L 2 xL 2 



= I L Jr™ /(* ~ s ' x ~ y)9(s,y)dyds\ L2y , L2 

(14) 

< \$ a a +1 !^\r s >- y (t,x)\\g{s,y)\dyds\ L2xL2 

< l/lL2 xL 2|g| L i xL i<is 

where f~ l ~ y {x) '■= f(x — y) denotes the function / shifted by —y, and we 
may use Minkowski's inequality. The functions considered are local with 
respect to time but we may consider s — > \/(s,.)\lp, t — > \g~ s (t, .)\ L i as 
functions in L 1 (M) by defining them to be zero in 1R \ [a, a + 1] . 

Now this type of Young inequalities may not be applied immediately in 
our situation, where one part of the convolution is a Gaussian. The appli- 
cation of Fourier transforms of fundamental solutions of the heat equation 
with respect to space and time variables shows this. For the partial (not 
normed) Fourier transformation with respect to the spatial variables we get 
for t > s 

j Rn exp(2^)| 1 exp (-j^) \dz 

(15) 

= exp (-4e(f - s)ir 2 ^ 2 ) . 
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Now as s 1 1 this equation becomes 1 (the formal Fourier transform of the 5- 
distribution) , and this is not in L . However for the truncated fundamental 
solution 



Mz)GS - s, z) := Mz) {2V J {t _ s))n exp ("^y) (16) 

the situation is different. Here 4>i 6 C°° (i?i(0))), i.e., with support in -Bi(O), 
and with 4>\{x) = 1 for \x\ < 0.5. Note that 4>i and l — (f>i build a partition of 
unity on M™. So the idea for estimating the increments 6v*' p ' ' * is to split up 
the integral of their representation and estimate one convolution summand 
with factor 

Gf(t-s,z):=Mz)Ge(t-s,z) (17) 
via the Young inequality and the other convolution summand with factor 

Gi l ~ B \t-s,z) := (l-Mz))G e (t-s,z) (18) 

via a weighted convolution estimates for L 2 functions. Indeed for the trun- 
cated Gaussian 4>i(z)G e {t — s, z) we have 

\Mx-y)Ge(t-s,x-y)\ = \<h(x- V) (2y jj {t _ s)r exp (~fe=§) I 

= \ M * -y)(t- sT-i 2 ^ (fe^) m ^ exp (-^ 

< \ M x - y)(t - (^) m j^j^ 6X P (-fe5) I 

(19) 

where 

1 1 / z 2 \ m ( z 2 \ 

C: =T i (5F(v^(t) cx H-ssJi >0 (20) 

is a constant depending on e, but finite for each e > 0. We may use this 
estimate locally for m = 1, i.e. we may use the upper bound 

\Mx-y)G e (t-s,x- y )\ < \M* - y)C(t - sr^j—t—^l (21) 

which is L 1 for dimension n = 3, because of the localisation and may 
be used with mixed Young inequalities. We start with L 2 estimates for 
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8v*' p ' 1 ' (r, .) = 5v*' p ' l ' k (T, .) where we define 
5v* ,p,l,k (r,x) =: 5v*' p ' l ' h,1 (T,x) + 5v*' p,l ' k,2 (r, x) + 5v* ,p,l,k,3 (r,x) 

■= Pi / R » f- E"=i ^i' P ''' fc ~ 1 ^ i &- — ) (s,v)Gi(t -s;x- y)dyds 

f- E^i^'^" 1 ^^- ) (s,y)Gi(T-s;x-y)dyds 

+Pl Ii-i Jk« /r» (£- R n( z ~ vj) x 



aft*'-''- 1 ftotfM-i , fl p y.'.fc-» < 

(22) 

We have a convolution with respect to time and with respect to the spatial 
variables. We do some estimates for the L 2 x L 2 -norm, or, for the squared 
|. \l2 xL 2 -norm. We may consider r — >• <5u*j' p ' '^(t, x) = lp-info'*/'' ' ^(t, a?) 
as an L 1 function on R where l[j-i,z] is the function which equals 1 on the 
interval [I — 1,1] and is on R \ [I — 1,1]. We implicitly assume this for 
convenience without changing the symbol of the function. We now apply a 
Young inequality to the first term on the right side of (|22p . We obtain 

I r *,p,l,kA f \ 1 2 



= M-iJr« (- EjU <foj' p ''** — ) {s,y)Gf{T-s;.-y)dyd 

^ W-iJk«I (E"=i*«J' / ''' , *~ 1 ^bJ — J (s,y)||Gf(r- s;. 

< IpA-i/ili/mn I (Ei=il* u J ,Mfc-1 l) (s,y)||Gf( r -s;--y)|^li2 



L 2 xL 2 



L 2 xL 2 



< I f^C fe _inmax ;7e{1) ... )ri} /^ 1 / Rn f |<5?;*' p ''' fc ^(r - s, x - y)\ \\Gf (s,y)\dyds) |^2 xi 2 



|2 

L 2 xL 2 



< |^Cfc_inmax je{1) ... in} \5v*' p " ' \ s > y (t,x) J R j Rn \Gf{s,y)\dyds\ 

< pfCl_ in 2 max je{1 ... in} \^)' p ' l ' k ~ l \ 2 L2xL2 \Gf\ 2 LlxLl 

< p 2 C 2 _ 1 n 2 (^) 2 max j6{li ... in} \6v^% xVi 

(23) 
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where 

Ck-i= sup (l + |y| 2 ) \D>*' P ' l ' k ~\s,y) , (24) 

se[l-l,l],yeR3 0<\a\<2 

where denotes the multivariate partial derivative of order a with multi- 
index a = («i, a2, as), and where 

C G '■= \ G f Li X Li- ( 25 ) 

Note that at this point the factor (|24|) seems obscur and we have could set 
it to 1 at this point. However, later we shall need weighted convolutions 
in I? for the estimation of the Leray projection term and then we need 
this wheight. The justification of the finiteness of Ct-i in (|24p . i.e., the 
polynomial decay of order m = 2 will be justified in the next section. An 
analogous argument leads to 

I r *,p,i,fe,l|2 

\ 0V (l~B)i lL2 xL 2 : — 



= \Pl Si-i /r« (" E"=i 6vj' p ' l ' k — \ {s,y)G[ 1 B \t-s;.- y)dyds 

< Pz 2 Ct 1 n2(^- B ) 2 max je{li ... in} \5v*/^\r, .)\\ 2 ^ 

(26) 

where 

Cg* '-=\G\- B \w (27) 
Summing up our result for the first term on the right side of (|22p we have 

(using (C|) 2 + (C£~ B) ) 2 < (C§ + C ( c~ B) ) 2 ) for all 1 < i < n 



2 

L 2 xL 2 



I r *,p,Z,fe,l|2 , 



G 



(28) 



x max je{1> ... )n} \5v*' p " ' \t, .)\ l2 |Gj| l1xL1 



<PiC k _ x n C G max je{1> ... in y\dv j \ L 2 xL 2, 



where 



C G = CB + & G B . (29) 
The latter estimate holds for all 1 < i < n, hence we have 



I r *,p,l,k,l 1 2 , 

max i6{1 ,... in} IffV" Il^xL 2 < 



< Pt C k-i n2c G max je{i,-,n} \Sv*' p ' l ' k \t,.)\ 2 l2xL2 . 



(30) 
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For the second term on the right side of ()22|) . i.e., for 5v* ,p,l,k ' 2 we get the 
analogous estimate 

I r *,p,l,k,2 1 \ 1 2 , 

max je{li ... ;re} \bv^> (t,)| l2 < 

(31) 

Finally we look at the third term on the right side of (|22p which is a 
double convolution integrated over time. We do the estimate in the special 
case n = 3 this time. Again we split up the function 

5v*^ k > 3 (r, .) = *,3*W»(t. •) + ^* { f^(r, .), (32) 

corresponding to summands with a truncated Gaussian and its complement 
as above. Each of this summand is again split up into two summands one 
of which corresponds to a truncated Laplacian kernel <j)\K t i (more precisely: 
its ith partial derivative), and the complement (1 — <t>i)K j. We define 

5v^(r, .) = 5v^(r, .) + 5v*/ {i %(r, .) (33) 

and 

c *,p,l,k,3 ( s _ , *,p,l,k,3 / s , c *, P ,l,k,3 , x /o 4 x 

where the references will be made precise in the following estimations. We 
start with the third term Sv*^ 1 ^' 3 which is defined via 

5v*/^' k ' 3 (T, .) = pi f Rn J Rn (faiz - y)£-K n (z - y)) x 

x(E^5^(^ + *^))(.,»)x <») 

JT j Gf(r — s, x — z)dydz 

This means that we have a truncated kernel (f)\{.)K^{.) and a truncated 
Gaussian Gf . The double subscript B indicates that we have bounded 
support in both cases. Recall that we 

setSv BBi ( T ,-) = for t 6 R\ [Z — 1, Z] 
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in order to write the time integrals conveniently. We obtain 

\ Sv BBi k ' 3 \hxL2 = /k/r» \Sv*/j£ k ' 3 (T,x)\ 2 dxdT 
= Jffi Jr« \PI fa« /r» (<M Z - y)£- R n( Z - Vj) X 

X I Lmj=l ^ "dxj 1 "d^j ) J l S ^y) X 

Gf(r — s, x — z)dydz\ 2 dxdr 

^ fm Jr* \pi Jr» I /r» " y)£- R n(z - y)J x 

(36) 

x (Em J= i ll^fe— + ^fe— I J (s,y)\dy\x 

\Gf' z (t — s,x — z)\dz\ 2 dxdr 

z)\dydz\ 2 dxdr 

where we recall our introduction of a function 0i £ C°° (-Bi(O))), i.e., with 
support in Bi(0), and with 0i(x) = 1 for \x\ < 0.5. Note that 4>i and 
1 — 01 build a partition of unity on W l . Furthermore we used our inductive 
assumption 

r% *,p,l,k— 1 o *,p,l,k—2 

sup I — -p. 1 ~ \<C k - X (37) 

The estimate (|36p is still not suitable since the function (j)i(.)-J^-K n (.) is in 
L 1 and not in L 2 . For this reason, we have to shift the x-argument to the 
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dSv 



,p,l,k-l 



function — 4^ — . We observe that the right side of ([36]) can be written as 

I Ir fan \pl Jr« fan (fi z+2v (x - y)£-Kn Z+2y (x - y)) x 



, f)/i,}*' p > l ' k ~ 1 , 

^[T, n m , =1 \-%^—\2C k ^)(r-s,y)\x 



\G Bl (s, z) \dydz\ dxdr 

= | f M J R n | Pi / R « J" R „ Ul Z+2V (y)£-Kn Z+2V (y)) X 



(38) 



x ( E^=i l^£^|2C fc _i )(T-8,x-y)\x 



\G Bl (s, z) \dydz\ dxdr. 
Next the function 

y^^ z+2y (y)^-K- z+2 y(y) (39) 

is certainly in L 1 for n = 3, so we can apply the Young inequality again. 
Hence, from (|3"6"|) we get 

i«:' M (r,.)ii 2xi2 

^4C|_ 1 (C|) 2 C^ 30i max mj - 6{1 ... n} 1 



9x m lL 2 xL 2 ' 



where 



C8 = \Gf\ LlxLl (41) 

as before, and where we use the same observation on mixed forms of gener- 
alized Young inequalities as above. Furthermore we used that for n = 3 we 
have 

\^Bm^K n {.)\ L1 <C K . i4>1 (42) 

for a finite constant Ck z 4>\ ■ 

Next for the second summand of fi v *' p ' l ' k ' 3 i. e for the summand of the 
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form ^j3(i_l#L we have by an analogous argument 

I r *,P,l,k,3 ,2 
\ 0V B(l-B)i\L 2 xL 2 

< p 2 i^U\! Rn (^ 2+2 % - y)l-Kn z+2y {. - y)) x 

x ^m,j=i| fe m I J l^yjay| i2xi2 |G ; | L i xL i| 

Sp^iC^) 2 ! Jr» (^*(y)£^*(v)) >< 



(43) 



< pfe c Li( C G B) ) 2n2c \^ max i,me{i r .,n} 1 8 ^'k m (r, -)|i2 



For the other two summands of ft v *' p ' l,k ^ ; i.e., for the summand of the form 
^ v *f-B)Bi an< ^ we have to estimate kernels of the form 

(1 - <h)Ki (44) 



while the treatment of the convolution with the Gaussian maintains, i.e., 
from the argument above we get 

I r *,p,l,k,3 |2 i I it *,p,l,k,3 |2 
l°' ; (l-_B)BilL 2 xL 2 l°' ; (l-_B)(l-_B)ilL 2 xL 2 

< P 2 i^UGl\ J Rn |(i - - y)^*^. - y)|x 

l aA,;^ I ' fc - 1 l \- , ' / . , |2 



x ( Eny=i I k m I J ( T ^)^L 2 
< p^Cf^C^max^^...^ | / Kn |(1 - <p 1 y z+2 y(y)x (45) 



w 9 &--^+2y/ \ I I 9Sv*' p ' l ' k ,- y . , 2 



/5?4C|_ 1 C2n 2 max ijme{li ... in} | J Rn |(1 - <^) z + 2 %)x 



^Kn z+2y (y)\j^ (l + M 2 ) \ 9 -%^—\-\r, x )dy\ 



Now we do no£ have (1 — (pi)(.)-£^-K n (.) € L 1 , and for this reason we 
introduced the constant C^-i in the form (|24|) above. This means that we 
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can give another factor = , T m to the convolution, i.e., we have 

I r *,p,Z,fc,3 |2 _|_ I e *,p,l,k,3 |2 

I (1— _B)Bi Il 2 xL 2 + l 0U (l- J B)(l- J B)jlL 2 xL 2 

< ^^Li^G^ma^meii,-,"} I Jr« K 1 ~ 0i)~ z+2j/ (y) x (46) 

x&^^^lTT^I^^r'^^^li 2 

with the same constant C^-i- Now we may use the fact that a function u is 
in L 2 if s > §, and for all 

u{x):= j{l + \y\ 2 )- s l 2 v{x-y)w{y)dy (47) 

for functions v,w £ L 2 , and such that 

Ml 2 < C s \v\ L2 \w\ L 2. (48) 
for a constant C s > 0. Indeed, we have 

{l-^){y)^-K n {y)eL 2 , (49) 

for n > 3 where we may denote an L 2 -bound by K 3L 2 , and 

— i (r,y)eL 2 . (50) 



Hence, from (1461) we obtain 



dx r 



I r *,p,Z,fe,3|2 _ is, *,p,Z,fe,3 |2 , I r *,P,l,k,3 ,2 

Il 2 xL 2 ~ l 0U (l-B)BilL 2 xL 2 + l 0U (l-B)(l-B)j Il 2 xL 2 

<-• „2/)r<2 r<2 r»2 T ,2r<2 rr , Q „ i d <fa J *' p ' i '' 1 1 / n|2 

Pl^k-l^G^K-iL 271 °s max j,me{l,- ,n} | fe^ V r i-J| L 2 



(51) 



Summing up and recalling that we have first order derivatives on the right 
side for some summands we write a H 1 norm on the right side (which suffices 
for our purposes). We have 



\K P ' l %^<P^C 2 _ lC 2 G n 2 : 



*,p,l,k—l 1 2 



(52) 



max j6{l,--- ,n} \$ v j' 



We note that for the old scheme considered in [3] and [T] we would have an 
estimate with right side in L 2 which is otherwise the same as in (|52p up to a 
constant related to the fact that we have one source term less and that the 
Gaussian estimates for the local fundamental solutions with variable drift 
term and their adjoints produce different constants. 
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4.2 I? x H 1 estimates and L 2 x if 2 -estimates 

In addition to the L 2 -estimates we need estimates for the first order partial 
derivatives of the components of the value function. We start with estimates 
for norms which include spatial derivatives, i.e., we start with I? x H 1 - 
estimates. At each stage k of the construction we may differentiate under 
the integral and start with the pointwise valid expression 

=■ ^K P ' l ' k '\r,x) + ^-5vT^\r, X ) + ^-5v*^(r,x) 

■= Pi ij-i J R n (- £?=i Sv* j ' p,l ' k ~ 1 ^^. — \ (s, y)G u {T -s;x- y)dyds 

+PiSi-iSan (~12j=iVj' P ' hk ~ 1 ^ i ^ ) (s,y)Gij{T-s;x-y)dyds 

+Pl J R n / Rn (£-K n (z - y)) X 

(dt)ll*' p ' l ' k ~ 1 /a *,p,l,k-l a *,p,l,k-2\\ 
£™,;=l (^fe— + ^fe— )) i",y)G u (r -8,x- z)dydzd S , 

(53) 

where the subscript denotes partial derivatives with respect to the jth 
spatial variable (as usual in Einstein notation). From the representation in 
(|53p we observe that the argument for the I? estimates can be repeated, if we 
have a L 1 x L 1 -bound for the first order partial derivatives of the Gaussian 
Gi (first order partial derivatives with respect to the spatial variables). We 
have to refine the simple estimate in (|19p a bit, observing that for a £ (1.5,2) 
and n = 3 we have 

i c «(* -.»-»)i= i , v ..;,.„,.-^' "p i 

< \( t _ s )a-l-n/2 1 ( (x-y? \ a 1 / (x-y) 2 \ I 

(54) 

< Iff - ^Q-l-n/2 1 / Qr-y) 2 \ m i i / (a-i;) 2 \ I 
- I V 1 S i |x-y| 2 «-i ^ (i-s) ^ (4e)« (2^)™ eX P ^ 4e(t-s) ) I 

<|C7(t- S )-i-"/a i5 ^ s=r |. 

It follows that we have local integrability with respect to time, since a — 1 — 
n/2 € (—1,0) and in space since \ x _^i a -i is locally integrable in dimension 
n = 3 for 2a — 1 6 (2, 3). Note that we have the same constant C as in the 
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simple estimate (|19|) above. Hence, we have 

\M^-y)G et i(t-8,x-y)\ LlxLl < Cg 1 , (55) 
and it is clear that 

|(1 - fa)(x - y)G e ,i{t -a,x- y)\ LW < cg~ B)1 . (56) 



Hence we may apply the same arguments as for L 2 x -estimates where we 
have to : 
and get 



have to replace the constants Cq for the Gaussian by Cq = Cq 1 + Cq 1 B ^ 



\ ov i Il 2 x//i _ 2-y=l \\ 0V Bi Il^xH 1 \ 0V (1-B)i \l?y.m 

< (T. + l)pf«t 1 (^)V(l + (C % + K^?(l + C*))x (57) 

I r *,p,l,k — 1 1 2 

Note that the additional factor n + 1 takes account of the fact that we need 
to estimate n + 1 = 4 terms with the method of L 2 -estimates above. For 
L 2 x ^-estimates we use convolution rules and partial integration to get the 
following representation of the second order partial derivatives of thevalue 
function. We have 



& 



-dv*' p " (t,x) 



:= « Jz-i Jr« (~ E"=i <^* W,fe — ) (s,y)G ltTn (r -s;x- y)dyds 
+Pl Il T -i /r" (~ Ei=i ^ W,fc_1 ^ i &- ) (s,y)Gi,m{r -s;x- y)dyds 

+PI Jl-! f R n J R n (^"^(z - ?/)) X 

x ^mj=i I m &, 1 ^te J J ( S '2/) X 



■J 



x Gi m (T — s,x — z)dydzds. 

(58) 

Here it becomes clear why we included second derivatives in the definition 
of Ck-i- Proceeding as before we need to apply the product rule in order 
to expand the derivatives j of the value functions above. This gives an 
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additional factor 2 at Ck-i- Furthermore we have 1 + n + n 2 terms that we 
have to estimate. Hence, 

I r *,p,l,h\2 _ v^3 ( ic *,p,l,k,j |2 _i_ it *,p,i,fc,j |2 i 

\ 0V i \L 2 xH2 - 2^j=l y\ ov Bi \L 2 xH 2 + \ 0V (l-B)i\L 2 xH 2 ) 

< (n 2 + n + l)pf4C«_ 1 (C^) a » a (l + (Cfc,* + ^ 3il ) 2 (l + C s 2 ))x (59) 

x max j6{1 ... i?l } |<ky" (r, .)| L 2 xH 2 

4.3 Proof of preservation of polynomial decay of order m, 
and higher order estimates 

In the representation of the functions v*' p ' k, \ 1 < i < n as in (|12p we have 
three summands 

PiJi-iJw (-E]=iS v * j ' P ' l ' k ~ 1 ^ 1 d^ — J {s,y)G l (T-s;x-y)dyds, 
Pi Ii-i / a » f- Ej=i uj'"'''*" 1 ^^ J (s,y)Gi(T - s;x - y)dyds, 

pi Ji-i t / K " - y)) x 

x I Lmj=l dx m \ dxj + Sajj J J K s iy) x 
xGi(t — s, x — z)dydzds. 

(60) 

All these summands have products of functions and derivatives of functions 
of the form 5v*' p ' l ' k ~ l , v * ,p,l > k ~ 1 known from the previous iteration step. Sim- 
ilar is true for representations of the functions v* ,p,k ' 1 ^ 1 < i < n. Hence 
if these functions 5v* ,p,l,k 1 ,v*' p ' l ' k 1 , and derivatives of these functions say 
of order up to |a| < m have polynomial decay of order p then products 
have polynomial decay of order 2p. This order of polynomial decay may be 
weakened by the convolution with the Gaussian or by the convolution with 
the Laplacian kernel in the Leray projection term, but we may expect that 
for p > n the polynomila decay of order p may be preserved. In the next 
lemma we analyze this. 

Lemma 4.1. Let r G [I — 1, 1] for some time step number I > 1. Polynomial 
decay of order m>2 and existence of continuous derivatives up to order m is 
preserved by the local scheme v*' p ' ' , 1 < i < n if polynomial decay of order 
m is given for the functions v*' p ' l ' k (r, .) G C m riH m , and 5v*' p ' l > k (T, .) G C m n 
H m . Especially, v*' p ' l > k+1 (r, .) G C m D H m , and 5v*' p ' l ' k+1 (r, .) G C m D H m . 
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Proof. We show that for m > 1 and < \a\ < 2 we have 

\D^6v*' p ' l ' k+1 \ < j^-, if |x| > 1 (61) 



a v *,p,l,k+l 



if this holds for | D®5v* ,p,l ' k \ . Similarly for higher order derivatives D 
with |a| > 2 and some p\ > 0, where D" = D^D" 2 ■ ■ ■ D" n denotes the mul- 
tivariate partial derivative with respect to the multiindex a = (a±, • • • , a n ). 
We have to estimate convolutions with the Gaussian G e for some small e > 0. 
The expressions for Q~ v *'P l > k+1 > l ' k+l anc [ and their derivatives (cf. (|58p 

for second order spatial derivatives; higher order spatial derivatives can be 
represented similarly with first order spatial derivatives of the Gaussian) 
involve terms which are essentially of the form 



or 



l-l 



l-l 



h(y)G t (t - s,x - y)dyds, (62) 



h(y)G e j (t- s,x- y)dyds, (63) 



where h is some function which is a functional of D%v*'P' l ' k {s, .) and I>£ 5v*'P' l ' k { 
for < |/3| < m, and where the latter functions are in C m r)H m and such that 
the functions themselves and their derivatives up to order m are assumed 
to be of polynomial decay of order m (according to inductive assumption). 
Furthermore we have a Gaussian factor of the form G t (t — s, x — y) (defined 
analogously as Gi above), or first order partial derivatives of this factor. We 
split up the integral of the convolution into two parts where one part is the 
integral for \y\ < On this domain we observe that the Gaussian has 
polynomial decay of any order. Indeed we have for < \t — s\ < 1, m > ^ 
and \y\ < Jfl > 



2 ' 



\G € (t - s, x - y) \ = \ 1 exp 



= K*- 5 )" 1 n/2 ( x Jy)2m ( ( ( t _s) ) (2v4f)« eX P ( _ 4e(t-s)) I ^ 
< \C(t - S ) m ~ n / 2 (a ._^ )Sm | < i^JT, 

where C = cC, where with z = ( x ~ v ^ we define 

•Jt—S 



C := sup I . ; — , ■ (z 2 ) m exp f — — ) I > 



(65) 



is a finite constant and c > is another finite constant. Furthermore, for 
the first order partial derivatives of the Gaussian we have a similar estimate, 
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i.e., we have for 1 < i < n, and < \t — s\ < 1, m > and \y\ < ^ > 
\G,A* -..«-»)! = I , v „;,.,„„ exp (- ggg) i 

- I f t ^m-n/2-1 l»-y| A^) 2 N | m 1 CXD f (x ~ y)2 ^ I 

(66) 

< |C(t - ^-n^-l ^^ l < ^ _ s)t n-n/2-l , 

for some constant C > 0, and with a locally integrable time factor which 
becomes nonsingular for m > 3 in the case of dimension n = 3. On the 
complementary domain |y| > ^j- we need some properties of the integrand 
h. We observe for \x\ > and some constants C, C > 

J«-i/{ M >M} $P G e(t~ s,x -y)dyds 

< Hi f{\ y \>i$}te{\ x -y\<i } ik G ^ ~ s > x ~ y^y ds 

+ III J{ ]y] >M} H]x „ y]>1} $? G & ~ s ^~ V)dy ds (67) 

< III /{M>M }&{ |^|< 1} &|C(t - -)- 1/2 (^l^ 

J{|j/|>^}&{l^-J/l>i} - 

Note if the functions v*' p ' ' , 5v*' p ' l ' k have polynomial decay of order m then 
all terms in the source function h except that Leray projection term have 
polynomial decay of order 2m hence for m > n the preceding observation 
indicates that the polynomial decay might be preserved. Hence the proof 
reduces to the observation that the integrand h in the form it has in the 
representation of Sv k +1 and their derivatives is of polynomial decay of order 
larger than p+n. Let us look at arbitrary partial derivatives of some maximal 
order m which are assumed to be of opolnomial decay of order m inductively. 
Now from (|22p we get for each 1 < j < n and a = (a±, • • • , ctj, • • • , a n ) =: 
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(3 + lj ■= («i , • • • > Pj + 1) • • • j «n) the representation 

: = Pi Ii-i kn (- E"=i 5v j' P ' l ' k ~ 1 ^ i ir j — ) ^ ( s ' y)°iA T ~a;x- y)dyds 
+Pi L-i /r« (- E"=i ^ W ' fc_1 ^fe- J ^ (8,y)Gi d (T -s;x- y)dyds 

+ PI Il-l fan fan (£-Kn(z ~ ^ 



x ^L m j=i m~ ^ m dxj 1 "fej J J i S ' X 
xG| j(r — s, x — z)dydzds, 

(68) 

where the subscript $ denotes multivariate partial derivatives with respect 
to the multiindex j3. Now concerning the first two terms in (|68p . i.e., 

D°dv*' p ' ' ^ (r, x) and D^8v*' p ' l ' k,2 (T,x) the respective integrands 



^ "flay 



3=1 



E „,,M-i^L_| (70) 

3=1 J 

are sums of products of functions Dx5v*' p ' l ' k 1 and D2v*' p,l ' k 1 of order 
|7M/5| < to, which are assumed to be of polynomial decay of order m. 
Hence all the integrands except for the integrand related to the Leray pro- 
jection term are of polynomial decay of order 2m, i.e., the argument above 
concerning the estimate for polynomial decay of order m for convolutions 
with Gaussians shows that D < £5v*' p ' 1, *'^(t,x) and D"Sv*' p,l,k ' 2 (T,x) are in- 
deed of polynomial decay of order m for \a\ < m. It remains to check the 
polynomial decay of order m for the term D"5v*' p ' l ' k,3 (T,x) for \a\ < m. 
Note that for \y\ < ^ we can use the Gaussian polynomial decay as above. 
Hence it is sufficient to estimate integrals of the form 



Ki(y - z)g(z)dz (71) 

{\v\>^} 
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where g is an integrand of the form 

*,p,i,fc— 1 n *,p,l,k — 2\ \ 

(72) 




dXj 



3 



with |/3 1 < m — 1. Note that 5 is again a functional determined by sums 
of products of functions Dx5v*' p ' l ' k 1 and Dlv* ,p,l,k 1 of order I7I, < m 
such that 5 is of polynomial decay of order m by assumption, i.e., we have 

b(*)l < ^ f ° r \*\ > t- ( 73 ) 

We can use a similar argument as above and write 
K,i(v ~ z )9{z)dz\ 



\y\>"4) 

^li"{|,|>M} & { N <M}^(^-%(^l 

+ 1 /{ H >M} & { N> M} " *M*H (74) 



+ 



This shows that D^5v*' p ' l ' k,3 (T,x) for |q| < m and m > n these functions 
are all of polynomial decay of order m, too. □ 

Next we consider the higher order estimates. Especially we need the 
H 1 x H m estimate for some m and with H 1 with respect to time. Product 
rules for Sobolev spaces allow us to reduce H l x .ff m -estimates to L 2 x H' m ~ 2 - 
estimates. Next the estimates are considered in the case n = 3. First we 
observe that for all k > and m > | we have 

o p,l,k 

^— e L 2 x H m ~ 2 , (75) 

OT 

and more generally for m > | + 2p we have 

9P 



5rP 



£ L 2 x H m ~ 2 P. (76) 



In this context note that the local-time functions q 1 tP — are considered to 
be trivially extended to the whole time as mentioned above. Consider the 
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first equation of (J2|) and the case k = 1. We have the representation 
— — - Pi v 2^j=i Pi 2^j=i v j 



(77) 

Knowing that v?' ' G L 2 x il m for the first term on the right side of (j77|) 



we have 

^E^^ 2 x^ m " 2 - ( 78 ) 

Furthermore, for v? ,l,k G L 2 x we observe that for m > | concerning the 
convection term on the right side of (|77|) we have for all r £ [Z — 1, /] 



- Pl J2vf' k -\r,.)^—(T,.) G fl™" 1 (79) 
j=i 

since the factors satisfy 7^''' fc 1 (r, .) £ i? m and 8 ^ (r, .) G iJ m_1 , hence 
the prodcut rule 

77 

|/S|ff» < C.I/Ih.Mh* for /,ff£F,s>- (80) 

applies for m > I in case n = 3. From our construction we know that the 
right side of (|77|) is locally continuous with repect to time. Hence, we have 

-^E^ 1 ^-^ 2 ^^ 1 ( 81 ) 

Finally, concerning the Leray projection term in (|77p we observe again that 
locally, i.e., for x — y G -B r (0) for some ball B r (0) of radius r > we have 

K ti G L\ (82) 

hence with an appropriate partition of unity, for example with (pi defined 
above we have 

G L 1 , and (1 - <fo)K ti G L 2 (83) 

for the first order partial derivatives of the Laplacian kernel. For m > | 
in case n = 3 the product rule ([80]) can be applied to the (relevant part) 
of the integrand in the Leray projection term, and using the pointwise rule 
ab < ^a 2 + T^b 2 in addition we conclude that 

/ „ p,i,fc-i d v p,l ' k ^ 1 \ 

ELi \ fe— M^niT- 1 . (84) 
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Hence for all r € [I — 1, 1] 



pi Ej,m=i /r« (M-- y)mTi K n(- - v)) TJj, m =i 



vra / dv 



o,l,k-l dv"/*- 1 



(r, y)dy 



+Pi Ejm=i /r» ( (i - - - ») 



n / dv 



X Ei,m=l 



(r, <E # 



m— 1 



(85) 

applying the product rule and different appropriate types of Young's in- 
equality to both summands (cf. also part II of this investigation). Again 
continuity with respect to time leads to 



Pi 52j,m=l Jr» dxi ^n(- V) ; 



A Z^j,m=l I fix, & 



(86) 



Hence, we have 



*,p,l,k 



(h 



e L 2 x F m - 2 



(87) 



for all fc > 0. The next step is to show that we have a contraction for some 
pi > 0. We observe 



dSv, 



,p,l,k 



< 



85v* 



p,l,k 



St 



+pi L,-=i 



L 2 xH m~2 

*,p,l,k—l d5v*' p ' l ' k ~ 1 



dxj 



L 2 xH" 



L 2 xH r 



< Pl v YTj=l 



d 2 8v* 



,p,l,k 



dv 



dx 'j L 2 xW 
*,p,l,k—l dv*' p ' l ' k ~ 1 



dxj 



L 2 xH r 



+PI 



Il-l J*R« J*R n (fai^niz y)j x 



x Em,i=l 



,p,l,k-l 



dx,. 



9x, 



L 2 xH r > 



We can estimate the right side of (|88p and for m > | we have the upper 
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bound 

pivn maxjg.fi,... ,n} 



8v„ 



,p,l,k 



+piC 5/2 n maxjg^... >n y 
+ftC 5 /2nmax i6{1 ... iIl} 6v' 3 
+PiC 5/2 C K n 2 maxj- m6{lj ... >n} 



L 2 xHm -2 
p,l,k-l 





dSv *,f,,l,k-l 




L 2 xH m-l 


dXj 


L 2 xH ,n-l 



,p,l,k-l 





Q v *,P,l,k-l 




L 2 xH m-l 


dXj 


L 2 xH m-l 



gsv *,p,l,k-l 



dx r , 



L 2 xH m-l 



dv 



*,p,l,k — l o *,p,Z,fc — 2 



+ 



L 2 xH m-\ 



From d88D and ([89]) we get 

dSv*' pJ ' k 



Or 



L 2 xH ,n-2 

-piC 5 / 2 nmax ie{1 ... in } 



< pivnmaxj^^...^ 
,p,l,k-i 



8v„ 



,p,l,k 



L 2 xH m-2 



6v. 



L 2 xH m-2 
,p,l,k-l 



r *,p,l,k— 1 

Sv 3 



+piC 5/2 nmax je{1) ... jTl} 

,mg{l,- ,n} 



L 2 xH m-2 



*,p,l,k—l 



L 2 xH m-2 



L 2 xH m-2 



ggv *,p,l,k-l 



dx r , 



L 2 xH m-2 



,p,l,k — l 



L 2 xH m-2 



+ 



*,p,l,k-2 



L 2 xH m-2 



Now form previous estimates we have the upper bound 



*,p,l,k— 1 



+ 



for some constant C k -i > 0, hence 

d6v*- p ' l ' k 



,p,l,k-2 



L 2 xH m-2 



<2C k . 



dr 



L 2 xHm -2 
2 



< pzz/nmax je{lj ... in } 



,p,l,k 



+PiC5/2^n + n 2 ) maxjgj^... >n y(l + C k )C k . 



L 2 xH m-2 



L 2 xH m-2 



(89) 



(90) 



(91) 



(92) 



Furthermore, from our previous estimates and for some p® we have contrac- 
tion for the first term on the right side of (|93p . Hence for some pi > 
independent of k we get 



,p,l,k 



L 2 xH m-2 



< jmax je{1 ... in} (l + C fe )C fe _ifo- 



*,p,l,k— 1 



L 2 xH m-2 



(93) 
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Similar construction estimates can be obtained by analogous methods suc- 
cessively for higher mixed derivatives the functions 

d5v*' p ' l ' k 

D«—^- (94) 

(for multiindex a), then for 

—^2—' ( 95 ) 

and then successively for higher order mixed and higher order time deriva- 
tives. 



5 Global linear bound of the Leray projection term, 
local and global solutions 

First let us consider the inductive construction of local regular solutions on 
[I — 1, 1] xK 3 by the scheme above. At each time step / > 1 having constructed 
v*' p ' l ~ l (l — 1, .) € C m n H m at time step I - 1 (at I = 1 these are just the 
initial data hi), as a consequence of the argument above we a have a time- 
local pointwise limit v*' p '\t, .) = v*' p ' 1 ~ 1 (t, .) + Efcli 5v*' p,l ' k (T, .) S H 2 for 
all 1 < i < n, where for n = 3 we have H 2 C C a uniformly in r £ [I — 1, 1]. 
Furthermore the functions of this series are Holder continuous with respect 
to time. Note that the first order time derivative -j^-v * ,p,l,k {r, .) of these 
members exist in H m as well for natural m > | as we proved in the previous 
section. We observed that v* ,p,l ' k € H m x H 2m can be obtained inductively 
for each m and this leads to full local regularity of the limit function of the 
local scheme. If we plug in the approximating function v* ,p ' 1 ' (r, .) into the 
local incompressible Navier-Stokes equation in the Leray projection form, 
then we have 

r\ *,p,Z,fc <-\0 *,p,l,k , , f\ *,p,l,k 

—97 Pi u 2^=1 dx j + Pi 2^=1 v j d x j 

-PI E I m =l /*. (&*n(* ~ V)) Zlm=l ( a C fc9 C!' fe ) ^ V)*V> 



Z^m=lJ*R« Ej,m=l f^ft^ <9z m 1 ( r ; 2/)<% 



+pi e • m =i (^n(x - y)) e i m=1 ( ^r^cr ) <?> V)*V 

(96) 
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As k f oo the right side of ([96]) becomes 

En <■ *,p,l,kdv*' p ' l ' k ^ n *,p,l,k d5v*' p ' 1, 

3=1 6V 3 ~fe^ + « 2Zj=l V j dxj 

Pi Ii-i fa« fan \mi K n(z - V)) x 

(ai,*,P,i,fe-l / *,p,l,k-l Q *,p,l,k-2\\ 
£mj=i (^fe— + ^fe— ) ) («, y)G.(r - «, * - 

(97) 

and since lim^oo 5v* ,p,l,k (T, x) = and lim^oo — — = for all (t, x) G 
[Z — 1, /] x R ra pointwise by our local contraction result, and since v*' p ' l ' k and 

OSv* 

— — are uniformly bounded, we observe that the right side expressed in 
(f§7|) goes to zero pointwise such that 

lim v*'"' l ' k e C 1 ' 2 ([Z - 1, Z] X R n ) , l<i<n (98) 

satisfies the local Navier-Stokes equation pointwise, and in a classical sense. 
Finally in order to show that the Leray projection term is globally linearly 
bounded, it suffices to show that the squared function 

l^\v*/\l,.)\ 2 H2 (99) 

grows linearly with the time step number Z. Now we can set up a controlled 
scheme similar as in pQ and in [2], i.e., to consider a scheme 

*,r,p,Lk *,P,i,k i I finr\\ 

v- K = Vj + rj (100) 

for some functions rj which have a uniform upper bound with respect to some 
regular norm \.\ H m x jj2m. Especially this upper bound should be indpendent 
of the time step number I. We shall see that it is possible to simplify the 
control function considered in pQ and [2]. Observe that the contraction 
results allow us to estimate the higher correction terms 



oo 

v *,P,i_ v *,P,i,o = J2Sv*' p ' l ' k (101) 

k=l 



in the local functional series on [I — 1, Z] x M. n at time step Z > 1 in terms of 
the function 

v*' p ' lfl -v*' p ' l ~\l-l,.), (102) 

i.e., the growth behavior of the latter function in (|102p determines the growth 
behavior of the former function (|101|) with respect to the relevant norm by 
the related contraction result. Assume for a moment that we can prove that 



v * >P> i,o_ v *, P ,i-i {l _ 1 ^ ) ^^ (103) 
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and assume inductively that 

max \v*' p ' l -\l - 1, .)\ 2 H2 <C + (l- 1)C, (104) 
je{l,-,n} 

which holds for I = 1 and for max je .n ... n i \hj\ 2 H2 for an appropriate constant 
for sure. Then assuming C > 4 w.l.o.g. we have 

max \v*' p ' l ~ 1 {l - 1, .)\ H 2 <^fC + {l- 1)C < i(C + (Z - 1)C). (105) 

je{i,-,n} J 2 

Now we have for some m > 2 

|5'y* , ' 0,i ' fe |2. mxH2m < imax i6{lj ... in} |^*' p ''' fc_1 |^ mxH2m , (106) 

where we may choose pi ~ y while for an appropriate choice of the constant 
factor (which transforms ~ to =) we have 



Cfe-i < (y/C + C(l-l) + l)foi all fc. (107) 
Assuming that we can realize the bound 



*, P ,l,o _ *,P,l-in _ 1 ) < (108) 

~ 2cVi 



we can choose such that we have indeed a contraction estimate with 
contraction constant 2C ^ - We get 

i *-p-U 

max je{l ) — ,n} Fj ' \H m xH 2m 

< max ie{lj ... in} |u*' p '' _1 (/ - 1, .) + E^o^j' P,/ ' P U m x//2m 

< max i6{1) ... in} \vj ,p ' l ~ 1 (l - 1, .)lir™xfl*» + E^lo l fo *' p ''~ 1,fe lff™x-H-2™ 

<Vc+(i-i)c + ^.. 



(109) 



Hence, 



max ie{li ... jn} \v*' p ' l (l, .)\ 2 HmxH 2 



<C + IC. 

(11 

Hence, if we can realize the estimate 



^W,o_^W-i (/ _ 1)) |<_J_ ) (m) 
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then the scheme is global. This leads to the idea to extend the scheme for 
*,p,l,k £ Q a controlled scheme for 

r,*,p,l,k *,P,l,k . I ;-,-in\ 

Vj' "* = Vj"" +r j , (112) 

where 

rj = -(v* j > p > l >°-v*> p > l - 1 (l-l,.)y (113) 

In this case we have to ensure that we can repeat the contraction estimates 
above with this simplified control function. Note that each time step I > 1 
we first compute v*' p ' l '° as before via 

' ^ - m» EU + pi T.U V 7 P ' 1 ^ 1 - L -)^(« - 1. •) = 

' pi e ■ m= i /,» (AM x - v)) Um=i (^^ 9 4^-) (i - i,v)dv> 

V*>P> 1 >°(1-1,.) =V*'P' l ~ 1 (l-l,.). 

(114) 

We then define with rj as in (|1 13[) 

r,*,p,l,0 *,p,l,0 i / /-,-|r\ 

?V = V + rj. (115) 

If we plug in this control function into the local Navier-Stokes equation a 
time-step I > 1 we get for each 1 < i < n two additional terms F [y r '* ,p ' 1 , 
and G (r') where -F is a bilinear functional of v r '* ,p '' and r', and G depends 
only on the control function. At each time step I > 1 we have to solve for 

r-, r,*,p,l no r,*,p,l , „ r,*,p,Z — 1 

< pi E",m=i /r- - v)) E - m =i ( ^g'" 1 ^!'" 1 ) - Ly)^ 

+ i? (v r .*.ft { ,r|) + G(r / ) , 
v*'^(i-l,.) =v*'* , - 1 (l-l > .). 

(116) 
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Then the higher approximations y^*' p ' l ' k are defined via the Cauchy problems 

f e\ r,*,p,l.k <-,o r. * ,p,l,k ,, . n r,*,olfe-l 

, pi E- m =i - y)) El m =i ^'a,':' ) (« - Ly)^ 

+ i? ( v r,*,Mfc-i jf J) +G(r') , 
v r,*,p,l,k(l - 1, .) = v r '*'' , ''- 1 (i - 1, .)■ 

(117) 

TT £ J.1 1 ■ I 1 i- ST T.*,p,l.k T,*,0,l,k r,*,p,l,k— 1 

Hence, tor the higher order corrections 0Vj = — fj we 

get the equations 

gg v ^*'P' > n d 2 5v T j , *' p ' l ' k ^ n r,*,p,l,k—l/j i \w 

— ^ w E J= i — fe| — + pi Ej=i *y (* - 1. •) x 

qc r,*,p,i,fe-l -, r\ r,*,p,l,k — 1 

x (« - 1, -) + Pi E"=i ^p' V - (i - 1, •) 

= Pi E], m =i fa« \fii K n{x - y)) x 

(dSv r " p ' l ' k ~ 1 /a r,*,p,t,fe-l a r,*,p,t,fe-a \ \ 

fowMQ, - 1,.) = 0. 

(118) 

Note that the G term cancels for the latter equations of functional incre- 
ments, and we have only an additional linear operator F (in a more general 
context we have defined this F explicitly in [1]). Now for this scheme the 
contraction estimates above can be repeated with rather trivial modifications 
for these linear terms, and we can realize the estimate (jllip . 

This leads to an alternative proof to the scheme with a control function 
as discussed pQ and [2], and it is at the same time a simplification of the 
scheme in [3], which are all based on the Leray projection formulation 
in [3]. The scheme in pQ and [2] with a control function have the advantage 
that the time step size of the scheme may be kept constant and the H 2 norm 
has a global upper bound for all time nonetheless. Numerically it may have 
also better stability. 

This scheme has to be modified similarly 
1 kampen@wias-berlin . de , kampenSmathalgorithm . de. 
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